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Abstract 

In this work, we report about the determination of nonperturbative OPE pa- 
rameters from fits of continuum perturbation theory to the Landau gauge quark 
propagator. The propagators are computed numerically using lattice QCD with 
Nf = 2 dynamical Wilson twisted mass fermions. We use four different values of 
the lattice spacing ranging from a ~ 0.1 fm to a ~ 0.05 fm as well as several quark 
masses per lattice spacing. This allows us to obtain continuum results for the chiral 
condensate and the up/down quark mass extrapolated to the physical point. The 
main results are the average up/down quark mass m^^(2 GeV) = 3.0 (4) (2) MeV 

at the physical point and {■lp^p)^^{2 GeV) = -(299 (26) (29) MeV)^ in the chiral 
limit. We have also studied nonperturbative contaminations of our results at small 
values of the momenta, which are often interpreted as the contribution of the gluon 
condensate {A'^). We do see contributions from such terms, which are, however, not 
stable over the order in perturbation theory. 

PACS numbers: ll.15.Ha, 12.38. Ge, 12.38.Aw 
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1 Introduction 



Quantum Chromodynamics (QCD) describes the strong interaction part in our current 
standard model of elementary particle physics. Due to its particular properties, predic- 
tions for fundamental parameters of QCD require nonperturbative methods. The main 
tool in this context is lattice QCD, which allows predictions from first principles. 
Also perturbative calculations in QCD, where applicable, play an important role for our 
understanding of QCD. For instance, in the regime of large momenta two-point functions 
{0{p) 0{p')) of some operator 0{p) can be written in terms of an operator product ex- 
pansion (OPE). The OPE contains coefficients carrying the dependence on the momenta, 
which can be computed in perturbation theory as asymptotic series in powers of the strong 
coupling Os, multiplied with local operators and appropriate factors of the quark mass. 
The matrix elements of those local operators, like e.g. the chiral condensate, are of purely 
nonperturbative nature. 

A most natural combination of the two aforementioned notions is to compute the mo- 
mentum dependence of two point functions nonperturbatively and compare these for 
large momenta with the OPE prediction. Such a method not only allows us to compute 
estimates for fundamental parameters of QCD, such as quark masses and condensates, 
but in principle also to determine the strong coupling constant a^. 

The quantity we consider in this work is the Landau gauge quark propagator P{k) in 
momentum space. Its OPE has the following form 



The coefficient functions Cxik"^), X = 1, m, ^/^t/^, . . . carry the whole momentum de- 
pendence and can be computed in perturbative QCD for large fc^. The nonperturbative 
information is encoded in the quark mass and the condensate (s). 

Since we are going to compute P{k) using lattice QCD, we cannot perform the calculation 
at arbitrary values for the momenta, since we have to fulfil the inequality 



The first inequality ensures small lattice artifacts and the second one the applicability of 
perturbation theory. It is not a priori clear that such a window exists and one of the ques- 
tions we try to answer in this work is whether the values of the lattice spacing available 
from state of the art lattice QCD simulations are yet sufficient for such an investiga- 
tion. Hence, we will pay special attention to both, lattice artifacts and nonperturbative 
contaminations of our results. 

Lattice artifacts have actually previously been a complication for the applicability of the 
method we are going to apply: when Wilson fermions are used naively in momentum 
space the leading contribution to the OPE of the quark propagator is a constant term 
proportional to the lattice spacing a. Even though this term will vanish eventually in 
the continuum limit, at finite value of a it will dominate the OPE, since the mass and 
condensate contributions are suppressed by powers of l/k"^. This makes the determination 
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of the chiral condensate difficult, whereas the quark mass can be determined if several 
values of it are investigated [1]. This complication can be circumvented by using the OPE 
of the pseudo-scalar vertex [2] or by working in the x-space [3]. 

Within the lattice formulation we are going to use here - the so-called Wilson twisted 
mass formulation of lattice QCD [4] - such a term can be avoided due to automatic 0{a) 
improvement [5]. In addition, we are going to remove lattice artifacts of order g^a^ from 
our data as computed in one-loop lattice perturbation theory [6, 7]. 
The procedure outlined above is not the only way to determine the quark mass and 
the chiral condensate. They have been determined previously for instance from fits of 
chiral perturbation theory to the data for the pseudoscalar decay constant and mass, 
see for instance [8]. But the analysis we are going to apply represents an independent 
way to determine these important standard model parameters with different systematics 
compared to other methods. 

The main results of this paper are determinations of the average up/down quark mass at 
the physical point and of the quark condensate, both within the MS scheme at renormal- 

ization scale 2 GeV^ 

mf^{2 GeV) = 3.0 (4) (2) MeV (2) 

and having performed the chiral limit 

(^^/^)^(2 GeV) = -(299 (26) (29) MeV)^. (3) 

For the chiral condensate we quote here the value stemming from one of our fit strategies 
(fit B), as explained later. The other fit strategies give slightly different, but compatible 
results. We also discuss nonperturbative contaminations at small momenta, called two- 
dimensional gluon condensate (A^) and provide effective values for it depending on the 
order of perturbation theory taken into account in the fits. 

The investigation presented in this paper is based on gauge configurations as produced by 
the European Twisted Mass collaboration (ETMC) with Nf = 2 quark flavors of Wilson 
twisted mass fermions [8, 9, 10, 11]. We refer the reader to these references for all the 
details of the simulations. 

The paper is organized as follows. In Sec. 2 we discuss the quark propagator in pertur- 
bation theory and in Sec. 3 our lattice formulation. In Sec. 4 we present our analysis 
strategy and in Sec. 5 the corresponding results. We conclude with a summary. 

2 The Quark Propagator in Perturbation Theory 

In perturbation theory the quark propagator is known up to three loops [12]. Recently, 
the OPE of the renormalized momentum space quark propagator in Landau gauge has 
been performed in Ref. [13] in the MS scheme. The authors have included terms up to 
mass dimension three in their calculations. According to its Lorentz structure, P{k) can 

^Thc first errors arc purely statistical and the second ones systematic, respectively. 
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be written as 

p{k) = ^^s{e)i + ^^v{e). (4) 

Assuming that we can use the OPE, one gets the following expansion for the scalar and 
vector form factors S{k'^) and l^(fc^) [13] (up to operators of dimension three): 



+ p mg{A ) + (^^) 



(5) 



and 

= C,(k')l + ^^m? + ^^{A^) . (6) 



While the quark mass and the chiral condensate (ipip) are clearly of physical origin, the 
existence of a gluon condensate (A^) is debatable. However, in the OPE such a term is not 
excluded. The Wilson coefficients Ci, Cm, Ca^, Cm^, C^^, CmA^ and C^a are functions 
of the strong coupling constant ^^(/i). The expansion can be found in Ref. [13]. The 
value of as{fi) can be computed via renormalization group (RG) evolution as described 
in Ref. [14], which needs as an input. For this purpose we use a literature value, as 
discussed later on. 

Although known up to a certain order in a^, the perturbative series of the Wilson coeffi- 
cients in Eqs. (5) and (6) may be truncated further in order to study the systematic effect 
of the truncation. For this purpose we stop the evaluation of the perturbative series at 
the rimax'th power in including at maximum terms of order a"™'"': 

Cx = C^x + CWs + ••• + C'^r^-'a^-^ (7) 

where X stands for 1, A^, m^,.... The truncation is consistently done in the Wilson 
coefficients as well as in the evolution of as by means of the /3-function. 



3 Lattice Formulation 



The lattice quark propagator is calculated within the twisted mass formulation of Nf = 2 
QCD [8, 9, 10, 11]. For a review see Ref. [15]. The gauge action used to generate the 
ensembles was the tree level Symanzik improved gauge action (tlSym) [16], viz. 



Sn 



ii,i/=i 



6o 5^ {1 - ReTr (f/^^^^^ } + 6i {1 - ReTr (t/^-^y } 

— 1/12 and = 1 — 86i. The fermion 

^ ^ XxDtYaXx ■ (8) 



with the bare inverse gauge coupling (3 = G/qq, bi 
action in the so-called twisted basis is given by: 
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Here represents the lattice Wilson Dirac operator, ttiq is the usual bare quark mass 
and fig is the bare twisted quark mass, which is multiplied by the third Pauli matrix rs 
acting in flavor space. Twisted mass fermions are said to be at maximal twist if the bare 
untwisted quark mass mo is tuned to its critical value merit, the situation we are working 
in. At maximal twist, the twisted quark mass /x^ is related directly to the physical quark 
mass and renormalizes multiplicatively only. Many mixings under renormalization are 
expected to be simplified [5, 17]. And - most importantly - as was first shown in Ref. [5], 
physical observables are automatically 0{a) improved without the need to determine any 
operator-specific improvement coefficients. For details on tuning to maximal twist we 
refer the reader to Ref. [8]. The aforementioned twisted basis x^X is at maximal twist 
related to the standard physical basis ip, ip via the axial chiral rotation 

In this framework we compute the Landau gauge twisted quark propagator Ptm{x) in 
position space 

PUx) = {x,Xo) = {{D,^r')u, (10) 

where (...)(/ denotes the average over gauge field configurations {U} collected by the 
ETMC, which were gauge fixed to Landau gauge using the overrelaxation method de- 
scribed in [18].^ The quark propagators calculated on these configurations are Fourier 
transformed to momentum space and then rotated into the physical basis using Eq. (9) 
yielding 

P(/c) = ^ (1 + nsr^) Ptraik)^ (1 + nsr^) . (11) 

The details of the ETMC ensembles we used can be found in Table 1. In total we consider 
four values of the inverse gauge coupling /3 corresponding to values of the lattice spacing 
ranging from 0.1 fm to about 0.051 fm [8], with a statistics of 240 gauge configurations for 
most of the ensembles considered. The values of the (charged) pseudoscalar mass range 
from 600 MeV down to 250 MeV. For setting the scale we use the results published by 
ETMC in Ref. [8]. 

On the lattice, Eq. (4) is valid only up to lattice artifacts. In particular, since parity is not 
a good symmetry of Wilson twisted mass fermions at finite values of the lattice spacing, 
a parity violating term in the propagator is allowed, which is formally 0{a). Ignoring 
higher order lattice artifacts including 0(4) symmetry breaking terms, the lattice quark 
propagator can be written as follows [20] 

P{k) = -^^V{k') + ^^S{k') + z75r^^G(A;2) . (12) 

The parity violating term G{k'^) comes with opposite sign for up and down propagators. 
Averaging over up and down propagators hence eliminates these artifacts [20]. Note that 

^Gribov copy effects have been found to be small for large momenta in studies of the lattice gluon and 
ghost propagators in Ref. [19]. We have also carried out a study on the Gribov copy dependence of the 
quark propagator on a test ensemble with smaller lattice size and have found no ambiguities there. 
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Ensemble 




L^xT 




a [fm] 


A^conf 


^1 


3.80 


243 X 48 


0.0060 


^ 0.10 


231 








0.0080 




240 


A3 






0.0110 




240 








0.0165 




240 


Bi 


3.90 


243 X 48 


0.0040 


^ 0.085 


240 


B2 






0.0064 




240 


B3 






0.0085 




240 


B4 






0.0100 




240 


B, 






0.0150 




240 


Ci 


4.05 


32^ X 64 


0.0030 


^ 0.067 


240 


C2 






0.0060 




161 


C3 






0.0080 




163 




4.20 


48^ X 96 


0.0020 


^ 0.051 


188 






323 X 64 


0.0065 




200 



Table 1: Details of the ETMC gauge ensembles, used in the analysis. See Ref. [8] 
for more details. 



the parity violating 0{a) term encoded in G{k'^) contributes to S{k'^) for standard Wilson 
fermions. 

Rotational symmetry is broken at finite values of the lattice spacing. In order to reduce 
the impact of those artifacts on our results, we carried out the following steps: firstly, 
we use for the analysis the lattice tree level momenta of the quark propagator, = 
^sm{ak^ + with k^ = (27rn^)/(aL^). Secondly, we restrict ourselves to momenta 

lying near the lattice diagonal (i. e. off-axis) by applying the so-called cylinder cut [21] 

where iV^ = 0.5 -(1,1, 1, 1) is the lattice diagonal. The constant has been set to Ccyi = 1.6 
for the large lattice at /3 = 4.2 and to Ccyi = 1.1 elsewhere. Following this procedure, the 
data for S^k"^) and V{k^) show a sufficiently smooth behavior. 

Finally, we correct the form factors for lattice artifacts of order O^gQO^). These corrections 
have been computed explicitly in Ref. [6, 7] within one-loop lattice perturbation theory. 

4 Analysis 

The basic quantities considered in our analysis are the bare form factors V^k"^) and S{k^) 
of the lattice quark propagator Eq. (12) after applying the cylinder cut and after removing 
lattice artifacts of the order 0{gQa^). We shall denote these corrected form factors with S 
and V . The 0{gQa'^) artifacts have only modest effect on the ratio S/V (after performing 



6 




Figure 1: Effect of the correction for 0{gQa'^) artifacts on S/V (left) and on V 
(right) for /3 = 3.8. The effect becomes weaker for larger /3 as expected. 



the cylinder cut) while V alone receives substantial corrections. As an example we show 
in Fig. 1 uncorrected and corrected data for aS/V and V for /3 = 3.8 in the left and 
right panel, respectively. At this (smallest) /3-value one sees dramatic corrections to V. 
However, as visible in Fig. 2, the corrections decrease as the continuum limit is approached. 
Our goal is to determine the renormalized quark mass and the chiral condensate. Looking 
at the OPE's Eqs. (5) and (6) it seems appropriate to study S^k"^), since the quark mass 
can be determined from its leading quark mass dependence and the chiral condensate is 
the leading contribution in the chiral limit. However, if one wants to avoid the usage of 
renormalization constants, it appears to be useful to consider the scalar-to-vector form 
factor ratio 



(14) 



since the renormalization constant cancels out. S and V are then replaced by their pertur- 
bative series and the resulting expression is fitted to our numerical data for R{k'^) with the 
renormalized quark mass, the renormalized chiral condensate and possibly further terms 
as fit parameters. We remark in passing that as soon as the renormalization constants in 
R cancel in between nominator and denominator, the OPE of R is written in renormalized 
quantities only. Scheme and renormalization scale depend then on the scheme and scale 
the perturbative expansion is performed in. 

For the fits of the ratio R as well as of the formfactors S and V separately we have followed 
three different strategies. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 

(a) 



0.4 0.6 

(b) 



Figure 2: Same as in Fig. 1(b) but for /3 = 3.9 (left) and ^ = 4.05 (right). The 
X-axis always displays the same physical momentum range such that 
these figures may be compared. 



1. Fit A: 

The fits of Eq. (14) to the data are performed simultaneously to all ensembles at fixed 
/3- values, allowing for a different fit parameter value for each ensemble, but only 
for one global fit parameter corresponding to {ip'^^) ■ AH other contributions from the 
OPE are neglected. At each value of (3 the results for ruq are then interpolated to 
reference values of rpmps, tq denoting the so-called Sommer scale [22] and mpg the 
pseudo-scalar meson mass, respectively. Those interpolated results are extrapolated 
to the continuum limit and then to the physical point where appropriate. 

2. Fit B: 

Resigning to determine a value for the renormalized quark mass, we extrapolate 
our data for R{k^) to the chiral limit at each /3-value first. Only then we fit the 
perturbative series to the data at each /3-value, with terms proportional to powers of 
mq set to zero and still neglecting all other contributions from the OPE. Thereafter 
the result for {ipil)) is extrapolated to the continuum limit. 

3. Fits A' and B': 

At small values of k"^ our data is potentially contaminated by additional nonpertur- 
bative terms in the OPE. One example is what is often called the gluon condensate 
{A^). We shall avoid here the discussion of its mere existence (see also [23, 24, 25]), 
but we shall investigate whether our data is contaminated by effects that may ef- 
fectively look like the gluon condensate. Since including such a term into the fits 
of type A oi B applied to the ratio R appears to be not stable, we first determine 
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(A^) from t^(A;^) only, and use it afterwards as an input for a fit of type A and B. 
4. Fit C: 

As a further check of the sensitivity of the extraction of the quark and gluon conden- 
sates against different assumptions in the fitting procedures, we resign to determine 
the quark masses rUq (as in fits B and B') and instead we fix their values to the ones 
of the renormalized quark masses, which at maximal twist are given by the twisted 
bare quark masses, a fig, multiplied by 1/Zp, where Zp is the renormalization con- 
stant of the pseudoscalar current^. In what follows, we will refer to such values as 
the Zp-based quark masses. 

Contrary to the previous fits A and B, we do not use the data for the ratio given 
by Eq. (14), but instead the data for the two form factors, S and V, separately in 
order to improve the sensitivity to the value of the gluon condensate, and moreover 
we consider simultaneously all the data for the four values of the lattice spacing. To 
this end we have to consider explicitly the values of the renormalization constant 
of the quark field, Zg, which will be treated as free parameters, and we introduce 
also simple discretization terms, proportional to the square of the lattice spacing, 
for the quark masses and the condensates (see next section). 

We perform fully correlated fits using the inverse covariance matrix in our x^-functions 
as described in Ref. [27]. In order to estimate the systematic error induced by a specific 
choice of the fit range we have performed all fits in several fit ranges. The results we give 
are then weighted according to the x^-distribution function and our final result consists of 
the weighted average over all fit ranges [8]. The fit ranges are consistently chosen among 
different values of the lattice spacings in such a way that the physical momentum range 
in terms of Tq/c^ is kept approximately constant. The different sets of fit ranges we used 
are 

roVe [18,..., 58], [18,..., 64], [18,..., 66], 
[19,...,64],[20,...,64],[21,...,64], 

which in physical units lie within the range 

3.9 GeV^ <k^ < 14.6 GeV^ . (16) 

The statistical errors are estimated using a bootstrap procedure to propagate the errors 
consistently to the next step of the analysis. In order to estimate the systematic error 
related to the truncation in the perturbative series we carried out the fits using n^^x = 2 
and ninax = 3 and we take the difference in the results as systematic uncertainty. 
The evaluation of the perturbative series requires a value of as input, for which a value 
of is needed. For this purpose we use the value Aj^ = 0.330 (23) GeV [24] (see 
also Table 2) which is in good agreement with other Nf = 2 determinations of the same 
quantity [28, 29]. The error of this number is taken into account in our bootstrap analysis 
and contributes to the statistical errors of our fit results. 

■^The values of Zp in the MS scheme at a renormaUzation scale of 2 GeV can be read off from Ref. [20] 
at /3 = 3.80, 3.90, 4.05 and from Ref. [26] at /3 = 4.20. 
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m^o [MeV] rp [fm] [GeV] 
134.9766(6) 0.42(2) 0.330(23) 

Table 2: Physical quantities used as an input parameters for the fits and analysis. 

The values of m^o and tq have been taken from Refs. [30, 8]. The value 
of is taken from Ref. [24]. 
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Figure 3: Fits of R in lattice units to our data for ensemble Bi at f3 = 3.9 (left) 
and ensemble Di at 13 = 4.2 (right). The vertical lines indicate the fit 
range. These plots correspond to fit strategy A. 



5 Results 



5.1 Fit A: Determination of Quark Condensate and Mass 

Ignoring terms proportional to (A^) and higher order terms in the OPE Eqs. (5) and (6), 
we have performed fits to the chirally non-extrapolated data of S /V with the quark mass 
and the chiral condensate as fit parameters only. The fits are performed as discussed 
previously. Throughout this analysis the renormalization scale has been fixed at /i = 
2 GeV in the perturbative series. 

For all fit ranges quoted above the fits have produced stable results with acceptable 
X^/dof values ranging from 0.78 to 1.24. We could not increase the lower boundary 
beyond r^k"^ = 21 as the fit for /3 = 3.8 then turned unstable. In all considered fit ranges 
the fitted values of the quark mass and the chiral condensate have been compatible with 
each other within errors. 

In Figs. 3(a) and 3(b) we show exemplary fits for ensemble Bi at P = 3.9 and Di at 
/3 = 4.2 with dashed vertical lines indicating the chosen fit range. 
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In order to perform a continuum extrapolation of the fitted quark masses we first 
have to interpolate (extrapolate) the values to common reference points of the squared 
pseudoscalar mass mpg at each /3-value. We have chosen the following reference points 



rlmls e {0.49,0.81,1.21,1. 



60}, 



which allow us to use interpolations in the pseudoscalar mass in most of the cases. Only 
for P = 4.05 and /3 = 4.2 we have to perform a short extrapolation for the largest reference 
mass Tgrnpg = 1.60 and for /3 = 3.8 and /3 = 3.9 we have to rely on extrapolations for 
the smallest reference point. In Fig. 4 we show as an example the interpolation to these 
reference points for /3 = 3.9 and P = 4.2. 

For each reference point we perform a separate continuum limit of the quantity rorriq in 
which is shown in Fig. 5(a). The data appears to be compatible with a linear behavior 
in for all chosen reference points as expected. 

Finally, the continuum quark mass data has to be extrapolated to the physical pion mass 
m^o = 134.9766 (6) MeV [30], for which we use a linear curve with zero intercept (leaving 
the intercept free gives compatible results). The extrapolation is shown in Fig. 5(b). For 
(ipip) we have also performed a continuum extrapolation linear in as shown in Fig. 6. 
Note that we also tried to include chiral logs for the quark mass dependence of m|g, 
however, at our current precision this does not make a difference. 

After performing a weighted average over the chosen fit ranges we quote the following 
results for Fit A: 



where the first error is statistical and the second is a systematic error reflecting the 
uncertainties related to the fit range and to the truncation of the perturbative series. The 
systematic error due to the variation of the fit range is taken as the maximum deviation 
from the x^-^veraged result. This amounts to about 0.08 MeV ((5 MeV)^) systematic 
uncertainty for the quark mass (condensate). Furthermore, decreasing the perturbative 
order to n^^x = 2 results in a smaller fitted quark mass value as well as a higher value for 
(ipip)- The systematic errors we quote is the change of the central values when we apply 
this modification and added by the change caused by varying the fit range. 

5.2 Fit B: Determination of {i/jiIj) in the Chiral Limit 

As in the chiral limit the chiral condensate is the only nonperturbative parameter (dis- 
regarding again a possible gluon condensate) we expect it to be estimated more reliably 
and with less statistical error than in the finite mass case. Having calculated the quark 
propagator for two to five bare twisted quark mass values fig we can perform a chiral ex- 
trapolation for each lattice spacing separately. This limit has been performed linearly in 
in the bare twisted quark mass a/i^ for R{k^). Fig. 7(a) shows exemplary fits for (3 = 3.9 
at three representative values of a^k"^, one at the lower end, one in the middle and one at 



MS 



(2 GeV) 



(335 (37) (35) MeV)^ 



mf S(2 GeV) = 3.0 (4) (2) MeV , (17) 



11 




(a) 



PS 



(b) 



Figure 4: Interpolation to the four chosen pion mass reference points in m-pg for 
/3 = 3.9 (left) and (3 = 4.2 (right). We also show the linear fit used for 
interpolation to the reference points. The blue points correspond to the 
interpolated values. These plots correspond to strategy fit A. 



the upper end of the considered momentum range. For any lattice spacing and any other 
momentum not shown here the data is consistent with such an extrapolation. 
We have then fitted the mass extrapolated data via Eqs. (5) and (6) disregarding all other 
OPE terms. The corresponding one parameter fits in ('^'0) have produced x^/dof values 
in the range 0.8 to 1.8 and have been performed in the fit ranges: 

r2A;2G [7,...,64],[9,...,64],[12,...,64], 

[7,. ..,60], [7,. ..,58], [7,. ..,56]. ^ 

Note that we had to extend the fit range compared to fit A towards the infrared in order to 
be sensitive to the curvature of the chirally extrapolated data and to obtain stable results. 
As an example we show the fit for /3 = 4.05 in Fig. 8(a). The continuum extrapolation 
has again been performed in the lattice spacing squared and is shown in Fig. 8(b). After 
a weighted average over the different fit ranges and the continuum extrapolation we get 
the following result for the chiral condensate: 

^^;/ (2 GeV) = -(299 (26) (29) MeV)^ , (19) 

where again the second error is systematic. Lowering the perturbative order to rimax = 2 
results in a higher value of (ipip). The systematic error due to the use of different fit 
ranges is evaluated in the same manner as for fit A and amounts to (6 MeV)^. 
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0.02 0.04 0.06 0.4 0.8 1.2 1.6 

(a) (b) 



Figure 5: Continuum limit of rQijiq for the four chosen reference points (left) and 
extrapolation of the continuum extrapolated quark mass values rouiq 
to the physical pseudoscalar mass (right). The linear fit has been con- 
strained to go through the origin. This result again corresponds to fit 



5.3 Fits A' and B': Nonperturbative Contamination 

As discussed in the introduction, on the lattice we have to restrict the analysis of the 
quark propagator to a window in the squared momenta, as for too small momenta the 
perturbative expansion will not be valid and for too large momenta lattice artifacts will 
be too large. Of course, with any choice of this window one can never be sure to be free 
of this sort of artifacts. For this reason we extend our analysis by including the further 
terms m{A'^) and (A^) in the OPE. Unfortunately, it turns out that an inclusion of these 
terms as fit parameters in the fit to our data for R{k'^) appears to be not stable. This is 
why we follow a different strategy: we first determine an estimate for (A^) from V{k'^) in 
the chiral limit alone, and use this estimate to repeat the fits we discussed before. This 
is not a fully consistent treatment, but it should provide an estimate of the uncertainty 
in our results. 

In more detail, we first studied V{k'^) of the quark propagator according to Eq. (6) in 
the chiral limit. In this form factor the dimension two term represents the first non- 
perturbative OPE contribution when the quark mass is extrapolated to zero and thus 
one has a direct handle on this term. A similar study has been done using the inverse 
lattice quark propagator in Ref. [31]. Note that unlike R{k'^) the form factor V requires 
renormalization, which we include as free parameters into our fits. 

The chiral limit of V{k'^) has been performed constant in a^q. This is justified as can 
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Figure 6: Continuum limit of {ipil^) from fit A. 

be seen for instance in Fig. 7(b) where we show the chiral hmit of V for (3 = 3.90 for 
three values of a^/c^. We then fit the OPE formula, Eq. (6), to our data for V^(/c^) in the 
same fit ranges as used for the fits A and B. We obtain rather large values for x^/dof 
of about ^5 — 7 for /5 = 3.8 and (3 = 3.9 which are related to the fact that the cut 
in momenta we have applied does not work well for two points of V at a^k"^ ^1.3 and 
a^k"^ ^ 1.55 (see the slight spread of the data at these regions in Fig. 1(b)). If we had 
discarded these points we would have obtained x^/dof < 3.5 in both cases. The fits for 
the two smallest lattice spacings has yielded acceptable values of x^/dof staying below 
~ 1.3(1.5) for P = 4.05(4.2) (see Fig. 9(a)). Finally, the results have been extrapolated 
to the continuum limit linearly in a^, as shown in Fig. 9(b). For the gluon condensate in 
the continuum limit we obtain 

{A^)^ = 0.65 (09) (17) GeV' . 

The large systematic error is dominated by the truncation of the perturbative order. In 
order to better understand the role of the (A^) term in the OPE, we have studied the 
dependence of the values of (A"^) on the order of perturbation theory that has been used 
in Eq. (6). To this end we have truncated perturbation theory at the order a"™'''' and have 
performed fits with rimax ranging from 1 to 3. For these fits we have restricted ourselves 
to only one fit range with 9 < Tq/c^ < 64. The resulting values of {A'^) have then been 
extrapolated linearly in a^/rg to the continuum, as shown in Fig. 9(b) for rimax = 2. 
In Fig. 10 we show the continuum extrapolated expectation value of (A'^) as a function 
of nmax- From the figure we conclude that with increasing nmax the continuum value 
of {A"^) decreases with no saturation visible (yet). Hence, we might conclude that this 
dimension two term is effectively describing higher order terms in perturbation theory 
which are not included in our analysis, see Ref. [32] for a discussion. However, it seems 
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Figure 7: (left) Chiral limit for R sX (3 = 3.9. Data has been extrapolated assum- 
ing a linear dependence on the bare quark mass a^iq which is compatible 
with every momentum considered, (right) Chiral limit for 

V at P = 3.9 assuming a constant dependence on the bare quark mass 
a/ig. Here we follow fit strategy B. 



that this contamination is not negligible for our data when momenta as low as Tq/c^ ~ 9 
(A;^ ~ 2 GeV^) are included in the fit. 

Next we investigate the influence of this (A^) term on the results for (ipip) and the quark 
mass. To this end we use the values of (A^) determined as discussed above as an input 
for fit with strategies A and B with fit ranges given in Eq. (15) and Eq. (18), as explained 
previously. The result of fit B' is a slightly smaller value of (ipip) 

(2 GeV) = -(294 (25) MeV)^ , (20) 

which is however compatible within errors with the value from Fit B. Also the results of 
fit A' differ only slightly from fit A: 



Nt 



(2 GeV) = -(324 (37) MeV)^ , mf^{2 GeV) = 3.0 (4) MeV . (21) 



If we had performed fit A (no {A^) term) with the same fit range, we had obtained 
(2 GeV) = -(335 (37) MeV)^ , mf ^(2 GeV) = 3.0 (4) MeV . 
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Figure 8: (left) Fit of aK in the chiral limit using fit strategy B for /? = 4.05. The 
vertical lines indicate the fit range, (right) Continuum limit extrapola- 
tion of (^■0). 



Thus, the difference is well covered by the purely statistical error which is quoted here. It 
is worthwhile to note that in both cases we were able to get reasonable stable fits of the 
fitted parameters (^/'^/') and when including the (A^) term in the fit with fixed values. 
We conclude that although we were not able to fit both nonperturbative condensates 
independently we get consistent results for (V'?/') and rriq with and without the contribution 
of (^2). 



5.4 Fit C: Determination of the Quark and Gluon Condensates. 

A drawback of the previous fitting strategies is that the sensitivity of the scalar to the 
vector form factor ratio, Eq. (14), to the value of the gluon condensate is quite limited. 
As a matter of fact the gluon condensate is either assumed to be zero (as in the fits 
A and S) or fixed at the value extracted from the analysis of the data for the form 
factor V extrapolated to the chiral limit (as in the fits A and B'\ The main reason 
for such a limited sensitivity is that, by expanding the denominator in Eq. (14), the 
power corrections depending on (A^) appear always multiplied by the quark mass m^, 
which is a small quantity. Moreover, the results of the previous fits suggest that the 
values of quark and gluon condensates are anti-correlated to each other. Therefore, in 
fit C we use the data for the corrected form factors S and V ^ separately, in order to 
determine simultaneously both condensates. The price to be paid is the introduction of 
the renormalization constant of the quark field, Zg, which is treated as a free parameter 
for each value of the lattice spacing. 
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Figure 9: (left) Chirally extrapolated vector form factor V and fits of perturbation 
theory to data at /3 = 4.05. The vertical lines indicate the fit range, 
(right) The continuum extrapolation of rg(^^) for nmax = 2 in /tq. 



At the same time one of the main outcome of the previous analyses is that the discretiza- 
tion effects on the quark mass and on the condensates appear to be proportional to the 
square of the lattice spacing and, moreover, the light u/d quark mass obtained in the 
continuum limit and at the physical point turns out to be in good agreement with the 
existing estimates made by the ETMC (see Ref. [8]). Thus, we try also an alternative de- 
scription of the lattice artifacts by fixing the physical quark masses to their Zp-based 
values and by performing in Eqs. (5) and (6) the following replacements: 

rriq -)■ — (a/ig) ■ 1 + 

Zp a V ^0 

(V'^)^(^^)-(^l + D^^^^ , 

where D^, D^^ and D^-^ are free parameters (independent on the lattice spacing). Finally, 
two discretization terms of the form Ds(v)(^'^k'^ i where Ds{y) is a free parameter, are 
added to Eqs. (5) and (6), respectively, to take into account possible discretization effects 
proportional to the squared momentum. The impact of such terms on the extraction of 
the quark and gluon condensates turns out to be quite limited. 

A total of 640 data points are analysed using 11 free parameters, obtaining a per degree 
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Figure 10: nmax dependence of the fitted gluon condensate {A"^) after continuum 
extrapolation. 

of freedom of ~ 0.9. The quark and gluon condensates turn out to be: 



MS 



N 



f 

(^2)^13 ^ 0.56 (06) (12) GeV^ 

where the second error is the systematic one reflecting the uncertainty in the truncation 
of the perturbative series. The value obtained for (A^) agrees with the one used in fits 
A' and B'. Furthermore, our estimates for (A"^) agree with the result fi'^(^^)J^ioGeV ~ 
ara) GeV^ obtained in Ref. [31] from the analysis of the Landau gauge quark 



-(270 (15) (20) MeV)^ 



2 GeV, 



2.01 (11 

propagator as in the present paper. The latter, once evoluted at the scale fi 
corresponds to {A^)%GeV = 0-67 (04) (l":^^) GeV^ 
The result for the quark condensate as obtained from fit C is a bit below the results 
obtained in the previous fits. This means that different treatments of the quark mass de- 
pendence of the leading term of the OPE of the quark propagator may lead to a systematic 
effect of ~ 30 MeV on the extracted value of the quark condensate. 

Overall we conclude from the comparison of primed with non-primed fits that the con- 
tamination of nonperturbative effects in our data is inducing errors that are well covered 
by the uncertainties we quote. 

6 Conclusions 



We have presented a study of the quark propagator on the lattice using Nf = 2 Wilson 
twisted mass fermions. By comparing the numerical data to perturbative series we were 
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quantity final value fitting method 





-(335 (37) (35) MeV)^ 


Fit A 


ruq 


3.0 (4) (2) MeV 






-(299 (26) (29) MeV)^ 


Fit B 




-(324 (37) (34) MeV)^ 




rUq 


3.0 (4) (2) MeV 


Fit A' 




(0.65 (09) (17) GeV)2 






-(294 (25) (28) MeV)^ 


Fit B' 




(0.65 (09) (17) GeV)^ 




m/Nf 


-(270 (15) (20) MeVf 


Fit C 




(0.56 (06) (12) GeV)2 





Table 3: Final values of our fit parameters {ip^p), niq and {^), all in the MS 
scheme at scale /i = 2 GeV. The value niq corresponds to the physical 
point, whereas (V'V') and {A^) are understood to be defined in the chiral 
limit. We emphasize that (only) the fit strategies B and B' provide an 
explicite extrapolation of the data to the chiral limit. This is why the 
B-fit value for (V^V) is quoted in Eq. (3). 

able to determine estimates for the quark mass and the chiral condensate, which both 
are fundamental parameters of QCD. A summary of the results for the different analysis 
methods can be found in Table 3. From these results, we obtain our final estimates printed 
in bold font in Table 3 and quoted in Eqs. (2) and (3) and in the abstract. The results we 
obtain are well compatible with other lattice determinations obtained using alternative 
approaches, while the errors we quote appear to be larger. However, we believe that these 
errors give a fair estimate in particular of the systematic uncertainties involved in the 
kind of analysis we applied in this paper. 

Therefore, we conclude that a combined perturbative and lattice analysis of the quark 
propagator is possible with recent lattice data, even if the errors are still large. Smaller 
values of the lattice spacing are desirable, since they would allow us to include larger 
values of the momenta in the analysis. 

We have also studied nonperturbative contaminations of our results at small values of 
the momenta, which is in the literature often interpreted as a contribution of the gluon 
condensate A^. We do see contributions from such terms, which are, however, not stable 
over the order in perturbation theory. In fact, in the continuum limit the contribution 
decreases with increasing order in and we do not observe any saturation as is visible in 
Fig. 10. Still, our value for the gluon condensate is compatible with the findings reported 
in Ref. [31]. 
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